On Field-Induced Quantum Criticality in YbRh2Si2- 
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The quantum critical point (QCP) in YbRh2Si2 is an enigma for the itinerant view of QCP. In 
an alternative view, this QCP is intimately linked to the selective Mott localization of the heavy / 
electrons. Following a perusal of this unusual QCP, I study an Extended Periodic Anderson Model 
(EPAM) within DMFT. A quantum phase transition (FQPT), accompanied by a rapid change in 
the Fermi volume, is found near the quantum-critical end-point of the selective Mott transition in 
the /-electron sector. The theory accounts for a wide range of unusual, singular non-Fermi liquid 
features exhibited at this QCP in YbRh2Si2 in a natural way. 

PACS numbers: PACS numbers: 71.28-fd,71.30+h,72.10-d 



Quantum Criticality in solids and the potential to 
"tune in" to new, novel phases of matter in their prox- 
imity underpins a large cornponent of modern condensed 
matter physics research P, Q . The rare-earth compound 
YbRh2Si2 is a particularly enigmatic case in point, ex- 
hibiting a "strange" non-FL "phase" fanning out from 
a. T — QCP separating antiferromagnctic (AF) or- 
dered and heavy Fermi liquid (HFL) phases. YbRh2Si2 
is tuned to the QCP by minute magnetic field or chem- 
ical substitution. Experimental data indicate that the 
Hertz-Moriya-Millis (HMM) scenario Q does not account 
for the unusual responses in YbRh2Si2- An alternative 
view posits that these are associated with the destruc- 
tion of the Kondo effect itself, either by inter-site RKKY 
interactions 0, Q , or by selective localization of / elec- 
trons 0. Given that the "standard model" of /-band 
systems, the Periodic Anderson Model (PAM), empha- 
sizes the HFL aspect driven by quasilocal Kondo screen- 
ing 01 , these new observations call for mechanisms which 
destabilize Kondo singlet formation. In spite of vigorous 
attempts [H, Q , the problem is far from "being solved" . 

I begin by recapitulating salient features of the FQCP 
in YbRh2Si2'- 

(i) dc resistivity, p{T) ~ AT over three decades in T Q , 
(ii) specific heat, Cv{T) ~ T°-^ at very low T 0, (iii) 
anomalously slow (in frequency, uj) decay of optical con- 
ductivity, with linear- inw scattering rate, t~^{uj) ~ 
with a ~ 1 0, (iv) strongly T-dependent Hall constant, 
Rh{T), and cote h{T) ~ CiT^-^Ca, and (v) rapid change 
in the low-T value of Rh across the FQCP, extrapolat- 
ing to a jump as T ^ jlQj, suggesting a rapid change 
in the Fermi surface (FS) across the FQCP, (vi) static 
magnetic susceptibility, x(q = 0, T, i? = 0) ~ 7^-0.6 j^j. 
T > 0.3 K, along with large Korringa ratio indi- 
cating very strong ferromagnetic correlations close to the 
FQCP. And x{B) ~ {B - 5^)"°-^ scales with the A-co- 
efficient of the term in p{T) in the HFL regime, (vii) 
NMR derived Knight shift, Kg{T,B) and the relaxation 
rate, l/TiT scale with x(q 0, T, S = 0). 

(i) , (iii) and (iv) are reminiscent of what is seen in high- 



Tc cuprates in their "normal" state. Recently, (v) has 
also been seen in cuprates near optimal doping [l3| All 
these behaviors are at odds with the HMM theory 0], 
which predicts markedly different behavior 0]. So the 
elucidation of (i)-(vii) in a single theoretical picture re- 
mains a challenge. 

Here, I address these issues by proposing a modified 
PAM with extended /-hopping and hybridisation, as well 
as a direct coulomb interaction between the / electrons 
and conduction (c) electrons, dubbed Extended-PAM 
(EPAM). I study this EPAM using DMFT, showing how 
the non-FL behavior along a curve in parameter space 
is understood as a selective Mott localization, and dis- 
cuss how (i)-(vii) naturally follow therefrom. To the ex- 
tent that this non-FL behavior is tied to /-Mott physics, 
single-site DMFT should capture the relevant physics. I 
will also show how this non-FL state is unstable to either 
AF, or to a heavy FL (HFL) away from this curve, at 
T = 0. 

The Hamiltonian is H = Hq + Hi , with the band part 
described by 



<i,j>,CT <i,j>,cr 

and the local part, by 



<i,j>,a 



(1) 



(2) 

I take the c-band centered around E = and con- 
sider Uff = 00 (valid for / shells), so the / electrons 
are projected fermions, Xif^ = (1 — nif^„)fi„, sat- 
isfying [X.ic,X^.^,]+ = 5ij6„„'{l - Uif^^a). Using the 
Gutzwiller approximation, Xifu — Qafia with = 
(1 - n/)/(l - uf,) and n/. = (1/iV) E.(/l/..>. This 
implies (tf,ef) q^{tf,ef) and Vfc qaVfc in what fol- 
lows, and corresponds to the slave-boson mean-field the- 
ory (SB-MFT), yielding a narrow, coherent / band with 



2 



a width W 



K 



the mean-field Kondo scale. 



Consistent with LDA calculations [1J|, I take e/ 
Below, I investigate the fate of this SB-MFT Kondo scale 
in presence of strong, quantum fluctuations of the / oc- 
cupation, caused by the competition between mean-field 
coherence, {Vfc, {tf.p/Uff)) and incoherence, driven 
by C//C. 

I start by splitting the Vfc term 

Vtjt;) E<.,,>..(/lc,. + h.c) + vtjt;E 




h.c) and consider H = Hq + Hi with V^-^' = y/tftp to 
begin with. Using ajg = {ufia + vcia), big = (vf^^-ucia) 
with u = y^tf /{tf + tp), V — tpl {tf + tp). it is easy to 

see that H = Hq + Hi is Hq = -tY,< '""^ 



and Hi = Ufc J2i.a a' l^iacrn-ibcr' + £/ jltA'^iaa + n.iba + 



{a\^h„ + h.c)]. 

This is the spin S ^1/2 Falicov-Kimball model (FKM) 
with a local hybridisation term, which is finite whenever 
£/ 7^ 0. Remarkably, when = 0, this reduces to the 
pure 5 = 1/2 FKM! Below, we show how (e/ = 0, Vfc = 
y/tftp) separates two, different metallic phases. 

First, at — 0, we see that [nib,H] = for each 
i, implying a local U{1) invariance of H: local configura- 
tions with n;, = 0, 1 are rigorously degenerate. This is ex- 
actly the condition for having singular b- "number" fluctu- 
ations. As is known [15] , the symmetry unbroken metallic 
phase is consequently not a FL, but is dominated by a 
superposition of one-particle (rif, — 0) and two-particle 
{ni, — 1) states at low energy. The a-fermion propa- 
gator is, assuming a lorentzian unperturbed DOS with 
half-width W for analytical clarity, very simple, showing 
"upper" and "lower" Hubbard bands, with a pseudogap 
at Ef{^ 0): 



rib 



(3) 



The corresponding 6-fermion propagator has branch 
cut singular behavior precisely at Ep{— 0), lead- 
ing to singularities in the local, one- and two-particle 
responses: Pb(t^) — 9{uj)\uj\~^^~°'°'^ and Xab^ {^) — 
Jdte'^'{alb,,it);blaUO)) - e{Lo)\Lo\-(^'-'>~<l Here, 
a = {l/n)ta.n~^{Ufc/W) is the so-called s-wave phase 
shift of the Anderson-Nozieres-de Dominicis (AND) X- 
Ray Edge (XRE) problem. Obviously, the FL quasipar- 
ticle weight, Z = 0, and the w, T dependence of physical 
quantities will be governed by power-law responses. No- 
tice that, in the (/, c) basis, divergence of Xab" {oj) corre- 
sponds to extended (cf. non-local hybridisation), singular 
quantum fluctuations of the / occupation: it is precisely 
these fluctuations which destroy FLT at e/ = via the 
AND orthogonality catastrophe (OC) [lB|. Given that 
the 4/ 7 /2 level hybridizes with two "c" bands in real- 
ity [141, we get the OC exponent, a = 2ao. 

In our FKM with e/ 0, and at T 0, the total 
fermion number, n — Ua + Ub — rif + nc, jumps dis- 



continuously from n_ = (1/2) -I- {\/iT)tan~^{Ujc/2W) to 
n+ = (3/2) — {l/'K)tan~^{Ufc/2W) for a range of densi- 
ties, n, near unity. This corresponds to a sudden jump in 
the 6-occupation, rib, giving a first order "valence" tran- 
sition as £/ is tuned through /x(= 0). At finite T, this 
line of first-order transitions ends at a second order crit- 
ical end-point (CEP), and Ub varies very rapidly over an 
energy scale O^ksT) around e/ = 0, extrapolating to a 
jmnp T = 0. Notice that this jump in rif, depends on 
Ufc/t, which we choose henceforth to be such that this 
jump is vanishingly small p/7[ at T = 0, giving a quantum 
critical end-point (QCEP). This explicitly shows the link 
between emergence of singular non-FL behavior and the 
selective Mott localization of the 6-electrons. We empha- 
size that, with finite tf,Vfc, both, the f- and c-fermions 
remain mobile: only their combination, b„ = (vf^ — uCu), 
is localized. Recent slave boson approaches have pro- 
posed this "selective Mott transition" of /-electrons in 
the context of the QCPs in RE systems p]. Our work 
is a concrete, DMFT-bascd, realization of the "selective 
Mott" QPT, with a non-local hybridization. In contrast 
to earlier work [6], nFL behavior here arises from the 
AND-OC in the corresponding impurity problem as Vfc 
is varied across a critical value, V^^^ = ^/tftp. 

We now show how the unique observations at the 
FQCP in YbRh2Si2 are understood as a consequence of 
the AND-OC derived above. The singularity in the b- 
DOS implies that their contribution to thermodynamic 
responses dominates that of the "itinerant" a-fermions. 
Hence, the low-T specific heat is 



Ce((T) ~ T.limr,-,Qlm[GbbiuJ + "7)]U=5 



(4) 



giving the 7 co-efficient as 7(T) ~ T actually, 
goes like T^C^^^hogiT/Tcoh) as the QCP is approached. 
The log- factor comes from seeing that the DOS is ap- 
proximately a lorentzian with a maximum varying like 
T^(^^"' and a full width at half-maximum equal to 
airT. When Ep{~ 0) lies within this peak, we can 
write Ep{T) = —airT, whence the asymptotic form of 
Cei ^ r-(i-")log(T/rco/i) follows. The entropy is then 
directly obtained as S(T) = j{T')dT ~ T". With 
ao = 0.3, we thus find that both C(T), S(T) vary as 
with a — 0.6, in nice agreement with observations in the 
nFL regime as a function of T ^ . 

What about transport?. In the impurity 

limit (note that the lorentzian unperturbed DOS 
in DMFT will not modify the "impurity" re- 
sult), we have Gaoir) — {ttT pQ / sihIttTt)) and 
G'fco(r) = sgn{T)/2, whence the respective self- 
energies are ^^(t) Uj^Gao{T)GaoiT)Gao{T) and 
E6(r) = UjcGaoiT)GaoiT)GboiT). Direct evalua- 
tion followed by Fourier transformation then gives 

Sb(zw„) = ~liUfcPo)^[tOniHEp/T) - ^!{uJn/2TlT) - ttT] 

and Yia{i^n) — (1^^ + tt^T^). The dc resistivity within 



3 



DMFT is then pdc{T) ~ {m / ne'^)lmY.b{u})\^=T ^ AT, 
i.e, it is Hnear in T. The optical scattering rate will also 
be linear in lo. This is exactly in accord with observa- 
tions near the QCP in YbRh2Si2 [1, Interestingly, 
with log-singularities in Sf, above (which already imply 
Z = 0), higher-order terms, which must be carefully 
examined, lead precisely to the branch cut singular 
structure [13] for Pb{uj), characteristic of the "lattice 
X-ray edge" problem found in DMFT. 

To proceed, observe that the impurity model corre- 
sponding to HpKM can be bosonized in each radial direc- 
tion centered around the "impurity" site [20| . For general 
band-filling, n = na + nj, ^ 1 per site, the umklapp terms 
from Ufc are irrelevant and hence ignored. The bosonized 
Lagrangian then describes a collection of non-interacting 
charge- and spin density collective modes: 



^0-2^ — 



(r) 



-{drMr)f]dr 



(5) 

and Lo" = Z]p / dr(f>p{r)dr. Here, g,Up^„,Kp^a are 
explicit functions of Ufc/t. Thus, interactions simply 
"shift" the charge bosonic modes relative to their free 
values. Introducing the usual symmetric-antisymmetric 
(charge-spin) combinations of <j)p.cr{r), we see that the 
antisymmetric (spin) channel completely decouples from 
the charge channel: a kind of high-dimensional spin- 
charge separation! This has been strongly emphasized 
by Anderson |16j| in the cuprate context, and has impor- 
tant consequences, detailed below. 

6/ 7^ has two effects: (i) it moves the 6-fermion 
level away from Ep, and, (ii) finite a-b hybridisation 
generates a finite, but heavy 6-fermion mass, due to re- 
coil in the XRE problem [2l|, giving a small "coherence 
scale", erec = kBTcoh, below which HFL behavior ob- 
tains in the lattice model. The FL quasiparticle over- 
lap, Z ~ e-'='(*=°°), with C{t) = 2C/2jXiii^(i _ 
cos{ujt))dLLi. This gives Z ~ exp[U'j^{ln{K)/{l — k^))]. 
Here, n = ma/mi,, with the band mass of the a- 
fermion and the heavy mass of the b fermion. Hence 
Tcoh oc Z (note the difference from the SBMFT scale, 
T^^) increases with e/, as indeed observed in the re- 
gion to the right of the FQCP. In D — oo, the relevant 
hybridization, Cf J2i cr('^L^'o" + h.c) implies that the one- 
electron DOS will show a narrow, low-energy FL reso- 
nance, with upper/lower Hubbard bands at high ener- 
gies, as is known 15]. Away from Vfc — y/tftp, the 

term J2<i j> <j ^^fcifla^ja + h.c) also causes one-particle 
intersite hybridisation between the a, b fermions. This 
again gives the &-fermions a finite mass and results in 
another HFL with Tcoh oc Z << 1. Low-energy re- 
sponses are then those of a HFL with Z << 1: an en- 
hanced 7 = C,i{T)/T, x{T) = xo, Pdc{T) = po+AT\ etc, 
followed by a smooth crossover to the non-FL response 

tjt^, €f = 0. 



found for V^^ — 



Using the bosonized form, Eq.(7), above allows further 
progress in the nFL regime. Expressing the transverse 
spin correlation function as an average over the phase 
variables permits its evaluation using Lq^o-- The result, 
following [22], is 



A 



-F{-) 



(6) 



with Kp = y^vp/ivF + Ufc) and vp = "^t the Fermi ve- 
locity. F{x) is a scaling function, ~ x for a; << 1 and ~ 1 
for X >> 1. And — 1 for the SU(2) invariant case, 
but Kg^ < 1 including spin-orbit (s — a) coupling effects. 
This immediately yields the power-law T-dependence of 
the NMR relaxation rate as 



1 = -J2lmx+-{q,io) ~ 



(7) 



The uniform spin susceptibility follows as x(q = 0, T) ~ 
X-{K-^-K,)^ In DMFT, the singular-in-w part of x(q, uj) 
is independent of q. This explains why both T^^{T) 
and Knight shift, Ks{T), scale like r-(^^'-^p) like 
x(q = 0,T). With a reasonable choice of Ufc/t, we get 
Kp = 0.4, leading to very good agreement with the host 
of power-law behaviors found in the magnetic response 
near the FQCP in YbRh2Si2- In particular, with the 
choice = 1, x(q = 0,T),T^\T),KsiT) all follow 
a T^^-^ law. Further, taking ao = 0.3 (see above), we 
find x/l — T~'^ '^. Assuming b/T scaling, where b = 
(B—Bc) is the distance from the critical field, this implies 
x(&)/7(&) ~ and that xib) - b-°-^-f{b) ~ b-^-'^ 

near the FQCP. Further, with A{b) ~ 1/& [2], we find that 
the Woods-Saxon ratio, A/-j'^ ~ b~° '^ and A/x^ ^ b° '^, 
which is weakly 6-dependent and saturates at "higher" b. 
All these are experimentally seen 12|. While ao = 0.3 
as found in our model DMFT may change somewhat in 
a truly "first principles" theory, the qualitative theory- 
experiment agreement is compelling. 

The higher-Z? spin-charge separation implied by the 
bosonized form of the impurity model also leads to con- 
sistency with the magnetotransport results: the Hall re- 
laxation rate is now controlled by spinon-spinon scat- 



With 

T-\ 



tering, leading to cot0jy(T) ~ ciT^ + C2 [16|- 
Pdc{T) = AcT, we get the Hall resistivity, Pxy{T) 
Both are in good agreement with experiment. 

We now turn to the evolution of the FS across the 
FQCP in YbRh2Si2- Hall data suggest an abrupt recon- 
struction of the FS across the FQCP as T ^ 0. Thus, a 
large FS in the HFL regime, also seen in dHvA work [23|, 
abruptly goes over to a small FS on the AF side. Within 
DMFT, in the symmetry-unbroken metallic phase(s), the 
shape and size of the FS is not affected by interactions, 
since the self-energy is purely local: J^a.bik, u) — Ea,b(tL'). 
Exactly at the FQCP, i.e, at Vfc = y/tjt^^'^f = 0' ^^'^ FS 
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will be a single sheet with a volume corresponding to the 
a-fermion number. To the right of the FQCP, the finite 
SVfc, £/ gives a finite a — b hybridization, giving a HFL 
metal, as found above. In this regime, DMFT studies 
on the PAM with a relevant hybridization yield a large 
FS containing both, the lighter a- as well as the heavy 
&-fermions. Thus, the abrupt change in the FS volume is 
intimately linked with the selective Mott localization of 
the ba — {vfa — uc„) fermions; they "decouple" from the 
FS at the FQCP. This occurs exactly at the point where 
the FL coherence scale vanishes, giving a non-FL metal 
with low-energy singular responses. 

What about AF order? For small 6Vfc,e.f << 1, two- 
particle processes, generated to second order in 6Vfc, e/, 
are more relevant than the one-particle a — b hybridiza- 
tion. These processes couple two "impurities" , and lead 
to two-particle instabilities, as in coupled D = 1 Lut- 
tinger liquids To this order, extra terms, Hrll — 

-A'E<.,,>«l&..^l'a.v, with A2 ^ 0{{dVf,)yUf,), 
are generated in H. In D = oo, these are decoupled 
as H^il = -A2 E<^,J>JMabalbi^ + Mtala,^.^ + h.c). 

Solving H = Hp KM + H?l within DMFT should yield 
an AF metallic phase pj|: it will have the same sym- 
metry as the AF phase in an "iitinerant" view, since 
both have Ma — {a\^aj^a) > 0. However, I choose a 
different route. Bosonizing these terms in H^l/ , the sec- 
ond term, corresponding to AF order, generates a cosine 
term in the bosonized Lagrangian for the spin sector: 
^mt _ g^cos(/3</iCT), with P = y/SnKa- With spin-orbit 
interaction, Ka < 1 , and the Lagrangian in the spin sec- 
tor. 



La = ^o,dr + 91 J cos{(3(j)a)dr (8) 

is a quantum sine-Gordon model with a relevant cosine 
term. This leads to an AF ordered state, corresponding 
to a finite expectation value of the (jja field: {4>a) > [2^ . 
Thus, AF order here results as a particle-hole instabil- 
ity of the singular, non-FL metal derived above in the 
spin channel, rather than from a band FS instability, 
as would be the case in "conventional" cases where FS 
nesting features in a FL metal give itinerant magnetism: 
the latter picture cannot account for power-law responses 
seen at the QCP in YbRh2Si2- Thus, in our EPAM, at 
the FQCP, selective 6-fermion localization permits AF 
to arise simply due to "inter-impurity" (corresponding 
to onset of RKKY-like) b ~ b local moment correlations 
induced via "itinerant" a-fermions. Obviously, the FS 
now has a small volume, containing only the "itinerant" 
a-fermions. 

At V^^J, the AND-OC will always occur in the 
symmetry- unbroken metallic phase in any D < oo. Thus, 
we expect that our findings will survive inclusion of non- 
local correlations beyond DMFT. Also, the /-electrons 



are never strictly localized: only the 6-combination lo- 
calizes at vjl\ An ab-inito theory for a = 0.6, Kp = 0.4 
and Ka- used here is hard: here, we have employed plau- 
sible Ufc/t = 10 (this is the only free parameter in our 
model) values. A truly first-principles correlated pro- 
gram (e.g, LDA-I-DMFT) is required to derive them. We 
plan to address this issue in future. 

In conclusion, a local QCP, triggered by the AND- 
OC 24], is found in the DMFT solution of the EPAM as 
the model parameters are varied. Using high-D bosoniza- 
tion, non-FL responses with an uncanny resemblance to 
those found at the FQCP in YbRh2Si2 are uncovered. 
This QCP is unstable, either to a HFL, or to AF. All 
these findings are in very good qualitative agreement 
with the T — b phase diagram of YbRh2Si2, whose un- 
conventional QCP is thence proposed to be of the local 
type, and associated with the selective Mott localization 
in the EPAM. Our analysis is potentially applicable to 
other, d- and /-electron based systems showing non-FL 
behaviors near the T — > itinerant-localized transitions. 
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